arXiv:1504.05379vl [astro-ph.SR] 21 Apr 2015 


Astronomy & Astrophysics manuscript no. Martinez-Gomez,Soler,Terradas ©ESO 2015 

April 22, 2015 


Onset of the Kelvin-Helmholtz instability in partially ionized 

magnetic flux tubes 

D. Martmez-Gomez, R. Soler, and J. Terradas 

Departament de Fisica, Universitat de les Hies Balears, E-07122, Palma de Mallorca, Spain 
e-mail: david.martinez@uib.es 

April 22, 2015 


ABSTRACT 

Context. Recent observations of solar prominences show the presence of turbulent flows that may be caused by Kelvin-Helmholtz 
instabilites (KHI). However, the observed flow velocities are below the classical threshold for the onset of KHl in fully ionized 
plasmas. 

Aims. We investigate the effect of partial ionization on the onset of KHl in dense and cool cylindrical magnetic flux tubes surrounded 
by a hotter and lighter environment. 

Methods. The linearized governing equations of a partially ionized two-fluid plasma are used to describe the behavior of small- 
amplitude perturbations superimposed on a magnetic tube with longitudinal mass flow. A normal mode analysis is performed to obtain 
the dispersion relation for linear incompressible waves. We focus on the appearance of unstable solutions and study the dependence 
of their growth rates on various physical parameters. An analytical approximation of the KHI linear growth rate for slow flows and 
strong ion-neutral coupling is obtained. An application to solar prominence threads is given. 

Results. The presence of a neutral component in a plasma may contribute to the onset of the KHl even for sub-Alfvenic longitudinal 
shear flows. Collisions between ions and neutrals reduce the growth rates of the unstable perturbations but cannot completely suppress 
the instability. 

Conclusions. Turbulent flows in solar prominences with sub-Alfvenic flow velocities may be interpreted as consequences of KHl in 
partially ionized plasmas. 
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1. Introduction 


Recent observations of the solar atmosphere have shown the 
presence of turbulent flows in quie scent prominences (see 
Berger et'aDl2010l : lRvutova et al.ll2010ll . These phenomena have 
been interpreted in terms of the Rayleigh-Taylor instability (RTI) 
and the Kelvin-Helmhotz instability (KHI). The latter is a well- 
known hydrodynamic instability caused by a shear flow ve¬ 
locity at the interface between two fluids (see IChandrasekhaii 
196 ill . It is possible to find in the literature a great number of 
papers devoted to the study of this instabili ty in many astro- 
physi cal environments, such as Earth’s aurora (Hallinan & DavisI 
197(lll . protopla netary disks (iMichikoshi & Inutsukall2006h . the 
magnetopause dGuo et a n 1201 (Hi or planetary magnetospheres 
iOgilvie & Eitzenreiteiill989ll . In solar coronal plasmas this in- 
stability has been observed, e.g., in coronal mass ejections 
(lEoullon et al.ll2()Tl1l . 

Classical magn etohydrodynamic studies (see, e.g.. 


IChandrasekhad Il96lll have shown that, due to the effect of 
magnetic field, fully ionized incompressible plasmas are stable 
to small amplitude perturbations if the velocity of the shear 
flow is sub-Alfvenic. Accordingly, the magnetohydrodynamic 
KHI can only be triggered by super-Alfvenic shear flows. Some 
turbulent flows detected in quiescent prominences exhibit a 
behavior that resembles the non-linear stage of the KHI, but the 
measured velocitie s, lower than 30 km s * (IZirker et al.lll998l : 
Berger et al. 201 Oh. are below the threshold, >100 km s * (see 


Terradas et al.ll20C^ . to trigger this instability. So it would seem 


that those turbulences cannot be interpreted as consequences 
of KH instabilities. However, the condition mentioned in the 
previous lines only applies to fully ionized plasmas, and qui¬ 
escent prominences are not fully ionized but partially ionized, 
i.e., they are also composed of neutral particles that do not feel 
the magnetic field and therefore ignore its stabilizing effect. The 
existence of this neutral component may modify the criterion for 
the appearance of the KHI, allowing the onset of the instability 
even for sub-Alfvenic velocities. 


The KHI in par tially ionized incom pres sible plasmas has 
been studied by, e.g.. lWatson et alJ(l2004h an d lSoleretalJdloill 
and it has been found that neutrals are unstable even for sub- 
Alfvenic flows. Therefore, in the absence of certain stabilizing 
factors, like, e.g., surface tension, partially ionized incompress¬ 
ible plasmas are always unstable under the presence of a veloc¬ 
ity shear. But these results have been obtained for the case of 
a Cartesian interface. However, the magnetic field in the solar 
atmosphere is better represented by means of flux tubes. Proper¬ 
ties of waves in a fully ionized magnetic flu x tu be have been in¬ 
vestig ated by, e.g.. lEdwin & Roberts! (119831) and iGoossens et al.l 
(I 2 OO 9 I) . The effect of partial ionization in a cy lindrical geome¬ 
try has been studied bv ISoler et alJ (l2009l 20131). H owever, mas s 
flow was not included in the works bv ISoler et ^ (l2009ll20T^ . 
a necessary feature for the possible appearance of a KHI. Hence, 
the logical step forward is to include that previously ignored ef¬ 
fect. 


In the present work we use a multi-fluid theory to examine 
the influence of partial ionization on the onset of KHI in mag- 
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Fig. 1. Sketch of the model used in this work 


netic flux tubes cooler and denser than their enviro nment . More 
precisely, we use a two -fluid approximation (see iPriesd 119821: 
IZaaarashvili et al.l 201 ih that treats electrons and ions as a single 
fluid that interacts with the other component of the plasma, the 
neutral fluid, by means of collisions. We include a mass flow in 
the longitudinal direction of the flux tube that has a discontinuity 
at the interface that separates the two media of different densi¬ 
ties. In order to avoid further complexity in the model, we ignore 
effects like the surface tension of the fluids and compressibil¬ 
ity. Restricting ourselves to the linear regime, we superimpose 
small-amplitude perturbations to the equilibrium state and derive 
a dispersion relation for the incompressible waves generated by 
those perturbations. This dispersion relation is a g eneralization 
of the formulas found in lEdwin & RobertsI (Il983h for the fully 
ionized case. 

This paper is organized as follows: in Section 2 we describe 
our model and then present the basic governing equations. In 
Section 3 we derive the dispersion relation for linear incompress¬ 
ible waves and obtain an analytical approximation to the unsta¬ 
ble solution for slow, sub-Alfvenic flows and strong ion-neutral 
coupling. In Section 4 we perform a parametric study of the solu¬ 
tions to the dispersion relation and obtain the dependence of the 
KHI growth rate on the model parameters. In Section 5 we per¬ 
form an application of the theory to a solar prominence thread. 
Finally, the conclusions of this work are given in Section 6 . 


2. Model and equations 

2.1. Equilibrium state 

The equilibrium state is a partially ionized cylindrical magnetic 
flux tube of radius a embedded in an unbounded medium. We 
use cylindrical coordinates, namely r, ip and z for the radial, az¬ 
imuthal, and longitudinal coordinates, respectively. A sketch of 
the model can be found in Fig. [T] The subscripts '0' and 'ex' 
denote quantities related to the internal and external plasmas, 
respectively. The densities of ions and neutrals are pi and pn, re¬ 
spectively, and depend only on the radial direction as 


Pi(r) = 1 

Pi,0 

Pi,ex 

if r < a, 
if r > a. 

( 1 ) 

Pn(r) = < 

f Pn,0 

1 Pn,ex 

if r < a, 
if r > a. 

( 2 ) 

Hence, there is 

an abrupt jump 

in density between the internal 


and external plasmas. We have chosen the particular case when 
the internal plasma is denser than the external one. The mag¬ 
netic field, denoted by B, is constant and pointing along the flux 
tube axis, with the same value in both media, i.e., Bq - Bex- In 
addition, we consider a longitudinal mass flow with constant ve¬ 
locity denoted by U. The flow velocity is discontinuous at the 
boundary of the flux tube. 


2.2. Governing equations 

In this work we study the behavior of a partially ionized plasma 
using a two-fluid theory. We assume the plasma is composed 
of an ionized fluid made of ions and electrons, and a neutral 
fluid made of neutral particles. The two fluids interact by means 
of ion-neutral col lisions . The general two-fluid equations can 
be found in, e.g ., |Prie^ (119821) . IZaaarashvili et akf (1201 ll) and 
iKhomenko et al.l (l2014l) . Here, we restrict ourselves to the lin¬ 
earized version of those equations, which describe the evolu¬ 
tion of small-amplitude perturbations. Hence, the set of two- 
fluid equations that describe the behavior of linear incompress¬ 
ible perturbations superimposed on the equilibrium state are 

^ \ 1 

T + r i ^ - Vpie -H - (V X fl) X B - PnV„i (Vi - Vn) , (3) 

Ot OZj p 


dzj^" 


-Vpn-PnVniCVn-Vi), 


(4) 



b -V x(viXB), 


(5) 


V ■ Vi = V ■ v„ = 0. 


(6) 


In these equations, Vi and Vn are the velocities of ions and 
neutrals, respectively, pie and pn are the pressure perturbations 
of ion-electrons and neutrals, respectively, b is the magnetic field 
perturbation, p is the magnetic permeability, 7 is the adiabatic 
index and Vni is the neutral-ion collision frequency. In addition, 
we define the ionization fraction as, x — Pn I Pi- 

In the following calculations, we replace the velocities of 
ions and neutrals by their corresponding Fagrangian displace¬ 
ments, and given by 


dt dz 


(7) 


Vn = 


dt dz 


( 8 ) 


In addition, we define the total (thermal -H magnetic) pressure 
perturbation of the ionized fluid, P', as 


P' = Pie + 


B b 



(9) 


3. Normal mode analysis 

From here on, we follow the same procedure as in ISoler et ^ 
(I 2 OI 3 I) and perform a normal mode analysis. Since the equi¬ 
librium is uniform in the azimuthal and longitudinal directions, 
we express the perturbations as proportional to exp{iimp -r- ik-z), 
where m and are the azimuthal and longitudinal wavenumbers, 
respectively. We only retain the dependence of the perturbations 
on the radial direction. Furthermore, the temporal dependence is 
set as exp{-ia>t), where o) is the angular frequency. Combining 
equations ([3])-(|5]) we arrive at a system of four coupled equations 
for the radial components of the Fagrangian displacements, 
and and the pressures, P' and pn, namely 

dP' 

“T “ Pi ~ ^Pn^ni^^r,n5 (10) 
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dpn 

dr 




Pi (q^ - ojI) 


1 d(r^i-i) 
r dr 




+ i 


Vni 

Q. + iVni 



3.1. Dispersion relation and approximate KHI growth rate 

^2 2 ^ To find the dispersion relation that describes the behavior of the 
waves in this system we need to impose that P', p^, ^r,i and 
are continuous at r = a, i.e., at the boundary of the tube. After 
applying the boundary conditions we get a system of algebraic 
(12) equations for the constants Ai - A 4 . The non-trivial solution to 
the system provides us with the dispersion relation, namely 


/'A 2 2 \ 1 d(r^r,n) 


iVni 


Pn 


—r + kt 


Pn'^n 


Pi (q^ - tol') 


(Q + iVni)^ + iVni) 


Pn, 


(13) 


where Q. — 10 - Uk, is the Doppler-shifted frequency. Other pa¬ 
rameters that appear in the equations are the square of the mod¬ 
ified frequency, Q^, and the square of the Alfven frequency, a>\, 
defined as 




1 H- 


IXVni 
Q + iv„ 


and 


to' 


A = K<^A^ 


(14) 


(15) 


rjk-,a) K'^{k-a) 

:) 


^mik^Cl') Kii^ikvCi) (Oq -b i'Vni.o) (^ex "b Ib^ni.e 
X [b^ni.O (^ex "b Ib^ni.ex) “ b^ni.ex (^0 "b iv^i^o)] ~ 0 


(23) 


where the prime denotes the derivative of the modified Bessel 
function with respect to its argument. 

In order to simplify equation (|2^ . we use the so-called thin 
tube (TT) approximation, i.e., we assume k^a 1. We perform 
an asymptotic expansion of the modified Bessel functions for 
small arguments and m + Q and only keep the first term in the 
expansion. The resulting TT dispersion relation is 


where c\ is the square of the Alfven speed, computed as 


0 ^ 

= — ■ 


(16) 

ppi 



Now we combine the equations (fT0li-(fT3ll and 
coupled equations for the pressures, namely 

arrive at two un- 

d-P’ 1 dP' 1 

2 rii^\ , 

(17) 

Qf.2 f. Qf. ( 


d^Pn 1 dpn 
dr^ r dr 


(18) 


whose solutions are combinations of modified Bessel functions 
of the first and second kind, Imikp') and K^aik^r), respectively. 
We require the solutions to be regular at r = 0 and vanishing at 
r —> oo. Hence, 


P'{r) = 


Pn{r) = 



if r < a, 
if r > a. 

(19) 


if r < a. 

(20) 


if r > a. 


where A\ - A 4 are arbitrary constants. In turn, the radial com¬ 
ponents of the Lagrangian displacements of the two fluids are 
related to P' and pn as 


1 / dP' ^ . Vni dpn \ 

Pi _ (^2 ^ \ ' Q -t- iVni dr j ’ 


( 21 ) 


&,n = 


-I- l- 


PnD (Q -t- iVni) (D + iVni)^ ^f22 _ ^^2 ^ 

1 dP' 


dpn 

dr 


Q + iv. 


'>ni Pi (q2 - ojI) dr ■ 


( 22 ) 


[pi.O (ffo "b (yoVni.o) + Pi.ex ^Dex (Gex + (ttexVni.ex) e^)] 

^ [Pn.O^O (^0 + Ib'ni.o) "b Pn.ex^ex (^ex + Ib'ni.ex)] 

■b [Pn.O^Ob'ni.O + Pn,ex^exb^ni,ex] “ 0? (24) 


which is the same expression, w ith a slightly different notation, 
as Eq. (37) of ISoler et akl (12012h . obtained for the case of a Carte¬ 
sian interface. Hence, the geometrical effect associated to the 
cylindrical magnetic tube disappears in the TT limit. 

Let us consider the limit when the collision frequencies go 
to zero, i.e., when the two fluids are uncoupled and show a com¬ 
pletely independent behavior. If we neglect the terms associated 
with the ion-neutral collisions, we can recover from equation 
(l24li the dispersion relations for the classical hydro dynamic and 
magnetohydrodynamic KHI (IChandrasekhaii[l961h . In that case 
the dispersion relation for our model is given by 

[Pi,0(D^ - o) + - ^A^ex)] IPn.oD2 + Pn.exDLl = 0, (25) 

from which we can obtain separated solutions for the ionized 
component of the plasma and the neutral one. On the one hand, 
the solutions for the neutral fluid are given by 


u - k 


Pn,of7o -b Pn.exf^e; 
Pn,0 "b Pn.ex 


+ ik,\Uo - t/e: 


I '^jPnfiPn.&x 

Pn,0 + Pn.ex 


(26) 


where the branch with a positive imaginary part implies that the 
amplitudes of the perturbations grow with time. This growing so¬ 
lution exists for any value of the shear flow velocity. Therefore, 
the neutral fluid is always unstable in the presence of a shear 
flow. 

On the other hand, for the ionized component we have the 
following solutions: 


CJ - k 


+ k. 


Pi,of7o "b Pi.ex 

Pi,0 "b Pi,ex 

Bl + Bl 

p(pi,0 -b Pi,ex) 


(t/o-t/ex)^ 


Pi.OPi.f 


nl/2 


(Pi,0 + Pi,ex) 


(27) 
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from which it can be seen that the magnetohydrodynamic KHI 
only appears when 


|t/0-C/exl> 


^0 + ^ex Pi,0 + Pi,ex 


Pi,0Pi,ex 


A 


0 Pfex) 


Pi,0Pi,ex 


(28) 


The presence of this velocity threshold is caused by the effect 
of the longitudinal magnetic held on the ionized huid. The ve¬ 
locity threshold is super-Alfvenic, which indicates that the typi¬ 
cally observed hows along, e.g., prominence threads or spicules 
cannot trigger a KHI if the plasma is fully ionized. This same 
criterion for the velocity threshold can be found in other works 
devoted to the study of waves in fully ionized plasmas with lon¬ 
gitudinal mass hows. For instance, equa tion (l28Tl i s equi valent (if 

we assume CA.ex = 0) to equa tion (5) o fiRv utoval ([198^ _ 

In addition, according to iRvutoval (Il988h and iRvutova et al.1 
(1201 Ol) . a fully ionized plasma may be subject to an explosive 
instability if the shear How velocity. At/, fulhlls the following 
condition: 


ca.o 


p/2 

\ ATT 

- < AU < ca,o 

Pi,< 


Pi,0 Pi.ex 


1/2 


(29) 


Therefore, taking into account the suggested existence of the ex¬ 
plosive instability, three unstable branches would be expected to 
appear in the set of solutions to the dispersion relation: the hrst 
one, associated to the neutrals, would be present for any value 
of the shear How velocity; the second one, the explosive insta¬ 
bility, would appear between the shear how velocity thresholds 
given by equation ( l29b : and the third one, associated to the ions, 
would follow the criterion of equation (l28l l. (Note that if Bex ^ 0 
but Bex = Bq, as in our model, the lower and upper limits of 
equation (l29l) should be multiplied by the factor V2). 

Let us go back to the coupled case. The full dispersion re¬ 
lation, Eq. (l2Tt . must be solved numerically when ion-neutral 
collisions are at work. However, it is possible to hnd an approx¬ 
imate solution when the ion-neutral coupling is strong and sub- 
Alfvenic Hows are considered. For simplicity, we also assume 
that the external How velocity is zero, i.e., Lex = 0. For sub- 
Alfvenic Hows the only unstable solution we obtain from the 
dispersion relation in the uncoupled case is that associated with 
neutrals, Eq. (l26l l. since the magnetic Held is able to stabilize 
ions, Eq. (l27l i. Hence, we try to Hnd a correction to the neutrals- 
related solution due to ion-neutral collisions. To do so, we write 
to - tOQ + iy, where oiq is the neutrals’ unstable solution given 
in equation (l26l l and y is a small correction. We insert this ex¬ 
pression for ai into the TT dispersion relation, Eq. (l24li . and only 
keep up to Hrst-order terms in y and second-order terms in I/q. 
After some algebraic manipulations we Hnd a solution for y. We 
omit the details for the sake of simplicity. Finally, the following 
approximate solution for the frequency is obtained. 


kzUopnfi . 2^2 t/QPn,0Pn,ex 

co~ - +1-, -TT-T- (30) 

Pn,0 + Pn,ex (Pn,0 + Pn,ex){Vm,0Pn,0 + Vni.exPn.exl 


The approximated growth rate, i.e., the imaginary part of to, is 
quadratic in the flow velocity and inversely proportional to the 
collision frequencies. Hence, if the same parameters are consid¬ 
ered, the value of the growth rate is smaller in the strongly cou¬ 
pled case than in the uncoupled case (compare with Eq. (I26l 0. 
Therefore, ion-neutral collisions have a stabilizing effect. 


4. Exploring the parameter space 

Due to its complexity, equation (l23t must be solved numerically. 
In this section we study the dependence of the solutions of the 
dispersion relation with respect to various physical parameters. 
Also we compare the full results with the analytical approxima¬ 
tion shown in the previous section. 

Throughout this section we use dimensionless parameters. 
Unless otherwise stated we use p;,o/Pi,ex = 2, p„,o/Pn,cx = 2, 
ca.o = 1 and k^a = 0.1. For simplicity, we assume that the colli¬ 
sion frequency has the same value in both internal and external 
plasmas, so we drop the subscripts from Vni. In addition, we now 
focus on the kink mode, the only mode that causes displacements 
of the axis of the cylinder, so we use m - 1. All frequencies have 
been normalized with respect to the kink frequency, cuk, which 
is the frequency of t he kink wave in the TT limit in the fully ion¬ 
ized case (see, e.g.. lRvutov & RvutovalfT97^ISDruitlll98lh . The 
kink frequency is given by 




Pt0^A,0 

f P /,0 "I" Pi,ex 


(31) 


To begin with, we study the dependence of the solutions of the 
dispersion relation with respect to the shear flow velocity. Hence, 
we take AU = Uq- Lex as a free variable. We consider three dif¬ 
ferent values for the collision frequency, which allow us to inves¬ 
tigate the behavior of the solutions depending on the strength of 
the ion-neutral coupling. Fig. |2] displays the real and imaginary 
parts of the frequency as functions of the normalized shear flow 
velocity for a) weak coupling (Vni/wk = 0.1), b) intermediate 
coupling (Vni/wk = 1) and c) strong coupling (vni/wk =10). The 
red symbols represent the solutions obtained numerically from 
the complete dispersion relation, Eq. (|2^ . The full numerical 
results are compared with the analytical solutions in the strongly 
coupled limit (shown as blue solid lines) and in the uncoupled 
case (shown as blue dashed lines). In addition, the classical shear 
flow velocity threshold for the KHI in a fully ionized fluid (Eq. 
(I28l0 is denoted by the vertical dotted lines. The upper panels of 
Fig. 12 display the real part of the frequency. We observe a very 
similar behavior for the three studied cases. Initially, when the 
shear flow velocity is zero, we And two solutions with nonzero 
wr. These solutions are associated to the ionized fluid and corre- 
spond to the usual kink M HD wave found in fully ionized tubes 
(lEdwin & Robertsiri983l) . The solution with cur > 0 is the for¬ 
ward propagating kink wave, while the solution with «r < 0 
is the backward propagating kink wave. A third solution with 
wr > 0 emerges when the shear flow velocity increases from 
zero. The new solution is associated to the neutral component 
of the plasma, in the sense that this solution only appears in 
the presence of neutrals. However, note that such simple asso¬ 
ciations between solutions and fluids cannot be done when the 
coupling is high and ions and neutrals behave as a single fluid. 
As the flow velocity keeps increasing, the three solutions con¬ 
verge for a critical flow velocity that depends on the collision 
frequency. The stronger the ion-neutral coupling, the smaller the 
critical flow. From that point on, the real part of the frequency is 
proportional to AU and is well described by the real part of Eq. 
(l26l) or, equivalently, Eq. (l30l l. 

The lower panels of Fig.|2]display the imaginary part of the 
frequency. The differences between the panels are much more 
remarkable than before, meaning that the value of the collision 
frequency has an important effect on the imaginary part of the 
frequency. The shaded zone denotes the region of instability, i.e., 
toi > 0. In such a case, perturbations exponentially grow with 
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Fig. 2. Upper panels: tUR/tui; as a function of the normalized shear flow velocity, AUIca.o, for k.a = 0.1, m = 1, and three different collision 
frequencies (a) Vni/wi; = 0.1, b) Vni/tUk = 1 and c) Vni/<ttk = 10). Lower panels: as a function of AUIca,o for the same set of parameters 

as above. The red symbols represent the solutions of the full dispersion relation, i.e., Eq. iSJ; the blue solid lines correspond to the analytical 
approximation given by equation ( 130b and the blue dashed lines show the unstable branch of the neutral fluid when there is no coupling (Eq. Il26b ). 


time (in the opposite case perturbations are damped). We observe 
that for small shear flow velocities there is only one unstable so¬ 
lution, corresponding to that originally associated to the neutral 
component of the plasma. A second unstable branch (originally 
associated to ions) appears for larger flow velocities above the 
classical super-Alfvenic threshold (Eq. (l28l l). Importantly, the 
only instability present for slow, sub-Alfvenic speeds is that orig¬ 
inally associated to the neutral component of the plasma regard¬ 
less of the collisional coupling between ions and neutrals. Ion- 
neutral collisions reduce the growth rate of that instability to a 
great extent, b ut collisions are not able to completel y suppress 
the instability (I Watson et al.ll2004l ISoler et ^l2012h . We have 
overplotted with a blue solid line the analytical approximation 
to the growth rate (Eq. (l30l l). The approximation informs us that 
the growth rate is directly proportional to the square of the shear 
flow velocity and inversely proportional to the ion-neutral col¬ 
lision frequency. As expected, the approximation shows a good 
agreement with the numerical results for small shear flow veloc¬ 
ities. Eor the case of weak coupling, the approximation is rea¬ 
sonably good for flow velocities up to 40% the internal Alfven 
speed. When the collision frequency is increased, the range of 
agreement between the numerical results and the approximation 
is greatly extended to super-Alfvenic speeds. Einally, note that 
in the left panel there is a stable solution that is not present in the 
other two panel (see the bottom-most curve). The reason for this 
absence is that for high collision frequencies this solution moves 
out of the vertical scale used in the plots. We are not interested 
in this solution because it is always stable. 

Co ncerning the explosive instability hinted by iRvutoval 
(Il988h . for our set of parameters this phenomenon would appear 
in the range 2 < At//cA,o < Vb (when Vni = 0; the existence 
of collisions between ions and neutrals would modify this crite¬ 
rion). If we pay attention to the lower panels of Eig. we find 
that the only solution with positive imaginary part of the angular 
frequency is the one associated to the neutrals. Thus, we observe 


no evidence of the explosive instability. The only additional re¬ 
markable feature in that range is found in the upper panels of the 
same figure: the real part of the angular frequency of one of the 
solutions associated to ions turns from being negative to being 
positive, meaning a change in the direction of the propagation of 
the wave. 


To investigate in more detail the effect of ion-neutral colli¬ 
sions on the instability for slow flows, we perform the follow¬ 
ing study. We fix the shear flow velocity to AUIca,q = 1 and 
compute the frequencies as function of the ion-neutral collision 
frequency, Vni. The chosen flow velocity is bel ow the classical 
thresh old for the KHI in fully ionized plasmas dChandrasekhaj 
[MU), so that only the neutral component is unstable in this con¬ 
figuration. The results are displayed in Eigure[3l where the so¬ 
lutions originally associated to ions are shown as red diamonds 
and those originally associated to neutrals are plotted with blue 
crosses. There is always one unstable solution for any value of 
Vni, but its growth rate decreases when the collision frequency 
increases. The growth rate gets reduced because neutrals feel in¬ 
directly, through the collisions with ions, the stabilizing effect of 
the magnetic field. In addition, as discussed before, the analytical 
approximation for the growth rate, Eq. (l30l) . shows a good agree¬ 
ment with the numerical results for high values of the collision 
frequencies, as consistent with the assumptions behind the ap¬ 
proximation. Concerning the real part of the frequency, we note 
that the frequency of solutions asso ciated to ions decr eases until 
it reaches a plateau for Vni/wk > 1 dSoler et al.ll2013h . while the 
frequency of the solution associated to the neutrals stays constant 
all over the range. 

Eor the sake of completeness, we have also studied the be¬ 
havior of the solutions depending on m. The mode with m - 0 
is known as the sausage mode and produces expansions and 
contractions of the pla sma tube, without displacing its axis 
dEdwin & Roberts| 1983h. Modes w ith m > 1 are known as flut¬ 
ing modes ( Edwin & Robertsll983h . It must be noted that, as can 
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a) 




3 


q; 

3 



y / CJ 

ni ^ k 


5. Application to solar prominence threads 

In this section we perform a specific application to solar promi¬ 
nence plasmas. We solve the full dispersion relation, Eq. (l23t . 
with values representative of a quiescent prominence. Therefore, 
the internal medium represents a prominence thread with densi¬ 
ties of ions and neutrals such that pi o -i- Pn,o = kg m“^, 
a temperature of Tq = 7000 K and radius of a = 100 km; 
the external medium is composed of inter-thread plasma with 
Pi.ex + Pn.ex - 2x 10“'° kg m“^ and Tex = 35000 K, which cor¬ 
responds to the regime of prominence-corona transition region 
(PCTR). Densities and temperatures are chosen for the equilib¬ 
rium condition of the total pressure (thermal plus magnetic) to 
be fulfilled, i.e., the total pressure is the same in both media. The 
magnetic fields are Bq - B^x = 10 G. We focus on the kink 
mode, so we take m - The neutral-ion collision frequencies 
depend on the temperatures and densities and are computed us¬ 
ing the following expression (see iBraginskiilflOb^ 



Fig. 3. a) and b) oJilw'x for the kink mode (m=l) as a function 

of v,ii/tuit, with AUIca.o = 1 and k-a = 0.1. The red diamonds are the 
solutions originally associated to the ions when there is no coupling, 
while the blue crosses are the solutions for neutrals. The solid line is 
the analytical approximation given by Eq. OS). In b) the shaded area 
denotes the region of instability. 


be seen in Eq. (l24li . the TT limit is independent of the value of 
m for m 0; this fact implies that in the range of applicability 
of that approximation there will not be substantial variations in 
the behavior of the different modes. To observe some dissimilar¬ 
ities we need to choose parameters beyond the TT case. Thus, 
we keep the values of the parameters used in the previous sec¬ 
tion but change the dimensionless longitudinal wavenumber to a 
larger value, namely, k^a — 2, and take Vni/wk = 1- The results 
show that there are minor differences between each solution: the 
solutions weakly depend on m and, for large m, they become 
independent of this parameter. 

We have also repeated the calculations done in this section 
for higher values of the longitudinal wavenumber, k^. We have 
not included the results here because they are not significantly 
different from what we have already explained. When k~ is in¬ 
creased, the system distances from the regime of the TT approx¬ 
imation, but this only produces a slight variation in the real part 
of the normalized angular frequencies, while the imaginary part, 
i.e., the growth rates, remain almost unaltered. 


Pi liekBT 

" —A - 

Znip V mrip 


(32) 


where nip is the proton mass, k^ is the Boltzmann constant, 
(Tin ~ 5 X 10“'° m^ is the collisional cross section for a hydrogen 
plasma. 

Eig. ID displays the most unstable solution of Eq. (|2^ as a 
function of the shear flow velocity for three different values of 
the ionization fraction: the red dashed lines represent the fully 
ionized case (x - 0), the blue crosses represent a partially ion¬ 
ized situation Ctf = 4) and the black diamonds depict a weakly 
ionized case (x - 100). The left panel shows the results for 
a wavelength d = 100 km, that correspond to a longitudinal 
wavenumber k, = InjA - 2n x 10“^ m“'. In the right panel 
the wavelength used is d = 1000 km, so k^ -2nx 10“° m“'. The 
shaded zone of Eig. |4]denotes a range of typical velocities (from 
10 km s ' to 30 km s“') that have been measured in quiescent 
prominences (IZirker et alJI 19981: iBereer et al.ll2010[) . The limits 
of this zone could slightly vary depending on the observations 
that are chosen as reference, but this variation is not significant 
for our analysis. We see that for a fully ionized case the insta¬ 
bility only appears for shear flow velocities far from the detected 
values. On the contrary, the cases with a neutral component show 
instabilities for all the range of velocities. Hence, partial ion¬ 
ization may explain the occurrence of KHI in solar prominence 
plasmas even when the observed flows are below the classical 
threshold. 

If we compare both panels in Eig. 0] we notice that on the 
right panel the growth rates are lower than on the left one. More 
precisely, they are about one order of magnitude smaller when 
the shear flow velocities are large. Since the wavenumber on the 
right panel is one order of magnitude smaller than on the left one, 
this behavior is in agreement with Eq. (l26l l. On the other hand, 
when the velocities are low, the growth rates on the right are two 
orders of magnitude smaller, which is consistent with what it is 
derived from Eq. (l30l l. 

In addition, in Eig. |4]we have overplotted with solid lines the 
approximate analytical solutions given by the imaginary part of 
Eq. (l30l l. The growth rates obtained from the analytical approx¬ 
imation are within the same order of magnitude; thus, Eq. (l30l l 
may be used to calculate estimations of the growth rates of KHI 
in a real prominence thread without the need of solving the much 
more complex full dispersion relation. 

Our analysis have demonstrated that KHI may be present 
in quiescent prominences due to the effect of partial ionization. 
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Fig. 4. Application to solar prominence threads. Growth rates as functions of the shear flow velocity for the following set of parameters: Pi,Q +Pn,o = 
10“^ kg m^'^.pi ex +Pii,ex = 2X 10“'° kg m“^, Bq - Bex = 10“° T, fl = 100 km, Tq = 7000 K and = 35000 K; on the left panel the wavenumber is 
k- = InjA = Inx 10“° m“' and on the right panel k, = 2nx 10“® m“'. The red dashed lines correspond to a fully ionized plasma (y = 0), the blue 
crosses to a partially ionized case Cp = 4) and the black diamonds to a weakly ionized case (x = 100). The solid lines represent the solutions given 
by the analytical approximation in Eq. ( 130b . The shaded zone denotes the region of values of flow velocity that have been frequently measured in 
solar prominences. 


However, we are not yet allowed to state that this magnetohy¬ 
drodynamic instability may explain the observed turbulent flows. 
Before that it is necessary to check if the growth rates given by 
the theory are consistent with an instability that could be ac¬ 
tually observed. Estimated lifetim es of prominence threads are 
about 20 minutes (iLin et al.ll2005h . For the region of highest ob¬ 
served velocities we compute the following growth times of the 
instability (t s l/oji): a) A — 100 km ^ t ~ 100 s for the 
weakly ionized case (x - 100) and r ~ 1000 s for x = 4', 
b) A - 1000 km ^ t ~ 10"^ s (~ 2.8 hours) for^ = 100 and 
T ~ 10° s (~ 28 hours) for;p = 4. The growth times obtained 
for A = 1000 km are larger than the typical lifetime of a thread; 
therefore, an instability originated by a perturbation with this 
wavelength cannot be the cause of the observed turbulent flows. 
On the contrary, the growth times for the shortest wavelength 
are of the same order of magnitude or lower than the detected 
lifetimes. This means that during the life of a prominence thread 
there is enough time for the development of a KHI caused by a 
perturbation with that wavelength. 

In addition, the analytical approximation we have derived in 
Section [341 i.e., Eq. (l30ll. may be a use ful tool in the field of 
prominence seismology (lBallesteijl2014l) . We introduce a new 
parameter v, called mean collision frequency of the plasma, and 
defined through the following relation 

^ _ _ 2pn,0Pn,ex _ (33) 

k (Pn,0 + Pn,ex) (Pn.O^ni.O + Pn,exVni,ex) 

Thus, the imaginary part of Eq. (l30l l can be now written as 


Tkhi = 




V 42 ° 42 7 khi' 


(34) 


Values of the three parameters that appear in the right hand side 
of Eq. (l34l) . namely the flow velocity, the perturbation wave¬ 
length and the KHI growth rate, could be estimated from obser¬ 
vations. Consequently, through this formula, we could estimate 
the coupling degree between the two components of the plasma. 


6. Conclusions 

Although they may not embrace all the physics of the considered 
system, simple models like the one developed in this paper allow 
focusing on a particular effect and ease the interpretation of the 
results. 

Here, we have studied how the existence of a neutral com¬ 
ponent in a plasma affects the propagation of waves and the pos¬ 
sible occurrence of KHI in a cylindrical magnetic flux tube. We 
used a two-fluid theory to obtain the dispersion relation for lin¬ 
ear incompressible waves. Then, we studied the dependence of 
the solutions on several physical parameters, namely the shear 
flow velocity, the collision frequency between ions and neutrals, 
the longitudinal wavenumber and the azimuthal wavenumber. 
We have found that perturbations at an interface separating two 
partially ionized plasmas are unstable for any velocity shear, 
contrary to what occurs in fully ionized plasmas, in which the 
effect of the magnetic field only allows the onset of the KHI 
for super-AIfvenic shear flows. Perturbations with large longi¬ 
tudinal wavenumbers have higher growth rates than those with 
smaller values of that parameter. The two constituent fluids of 
the plasma, i.e., neutral and ionized, are coupled through col¬ 
lisions. This coupling holds a significant influence: increasing 
the ion-neutral collision frequency reduces the growth rates of 
the unstable perturbations, although it is not possible to avoid 
the onset of the KHI. These results are consistent with previ¬ 
ous works like the ones d evel oped in simpler co nfigurations by, 
e.g.. lWatson et al.l (l2004l) and lSoler et al.l (l20l^. Eg. of the 
present paper re duces to equa t ion (3 7) of ISoler et al.l (1201^ or 
equation (26) of lWatson et al.l (1200^ in the TT limit, i.e., when 
the product of the longitudinal wavenumber by the radius of the 
tube is much lower than one. Moreover, in the absence of an 
equilibrium flow and when the densities of neutrals go to zero, 
we can recover from o ur dispersion relation the incompressible 
limit of equation (8a) of lEdwin & Robertsl(ll983h . corresponding 
to a fully ionized cylindrical flux tube. 

We have also found that the solutions of the dispersion rela¬ 
tion slightly depend on the azimuthal wavenumber. Modes with 
a higher azimuthal wavenumber are more unstable, but this vari- 
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ation is only appreciable beyond the TT limit and it has not a 
great significance. 

Then, we applied our model to a thread of a quiescent promi¬ 
nence, by choosing values for the densities and temperatures typ¬ 
ical of those coronal features and solving the dispersion relation 
for several degrees of ionization and wavelengths. We have con¬ 
cluded that, for a certain combination of parameters, the turbu¬ 
lent flows detected in quiescent prominences may be interpreted 
as consequences of KHI in partially ionized plasmas. The growth 
rates of the instability increase when rising the ionization frac¬ 
tion, i.e., when the relative densities of neutrals increase. 

Furthermore, we have provided an analytical approximation 
of the KHI growth rates for slow shear flows and strong ion- 
neutral collisional coupling. This formula, Eq. (l30l l. is easier to 
handle than the full dispersion relation and thus easier to inter¬ 
pret; growth rates of the KHI show a quadratic dependence on 
the longitudinal wavenumber and the shear flow velocity, and are 
inversely proportional to the ion-neutral collision frequency. The 
analytical approximation may be useful in the field of promi¬ 
nence seismology. From it we can define a mean collision fre¬ 
quency, V (given by Eq. (l33T l). that provides an estimation of the 
coupling degree of the plasmas. Values of this parameter may be 
computed from observational data. 

And even though we have used a model where the magnetic 
field is equal in both media, the conclusions extracted from its 
analysis may also be valid when the magnetic field inside the 
flux tube differs from that outside. 

For the sake of simplicity, in our model we have ignored ef¬ 
fects like gravity, compressibility or surface tension. We have 
also chosen a particular alignment between the mass flow and 
the magnetic field: they are parallel, a configuration that gives a 
higher stability to the system. In future investigations the model 
used here can be improved including the mentioned effects, that 
may have an impact on the KHI. Moreover, we have focused 
only on the linear phase of the instability, while the non-linear 
regime would be also of interest. And we have assumed that the 
plasma is composed only of hydrogen, while actually more el¬ 
ements with several states of ionization are involved. A much 
more complex model that incorporates various of the mentioned 
refinements would be needed to fully understand the KHI in 
prominences, but such model could not be studied analytically 
and numerical simulations would be required. Our paper is an 
improvement from previous investigations but there is space for 
a lot of future work in this field. 
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